Mon. Not. R. Astron. Soc. 000. [THTTl C20041 Printed 2 February 2008 (MN WF&t style file v2.2) 



Inconsistency in theories of violent-relaxation 



I. Arad* and D. Lynden-Bell 

Institute of Astronomy, Madingley Road, Cambridge CBS OHA, UK 



in 
o 
o 

(N 

X> 
<D 
IX, 

(N 

m 
> 
oo 
(N 
t> 
OS 

o 
o 



Of 

6 

CO 

c3 



X 



Version of 2 February 2008 



ABSTRACT 

We examine an inconsistency in theories of violent-relaxation by Lynden-Bell and 
Nakamura. The inconsistency arises from the non-transitive nature of these theories: 
a system that under goes a violent-relaxation, relaxes and then upon an addition of 
energy, under goes violent-relaxation once again would settle in an equilibrium state 
that is different from the one that is predicted had the system would go directly from 
the initial to the final states. We conclude that a proper description of the violent- 
relaxation process cannot be achieved by equilibrium statistical mechanics approach, 
but instead a dynamical theory for the coarse-grained phase-space density is needed. 
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1 INTRODUCTION 

The main purpose of this paper is to demonstrate that the- 
ories predicting a definite statistical equilibrium, dependent 
only on the energy and the initial volume T(r))dr) at phase- 
space densities in the range 77 — ^ 77 + drj, do not predict the 
same final state when the system undergoes two violent re- 
laxation sessions separated in time, as they do when the two 
sessions are treated as one. This may be called an inconsis- 
tency or at best a lack of transitivity. 

Consider an iV-body gravitating system that starts in 
some state far from equilibrium, vibrates violently and set- 
tles to a dynamically steady state si that lasts long enough 
that it may be considered the final product of violent relax- 
ation with energy E\ . Now suppose that this system suffers 
a significant tidal disturbance from a passing object that 
causes violent vibrations and leaves out the system to relax 
again but now with energy E2- There are now two ways of 
predicting the outcome. Either we take the function to (77) 
giving the volume at each phase-space density from the ini- 
tial state so, or we use Ti{rj) the predicted outcome for that 
function after the first relaxation process. Of course if we 
used the fine grained phase-space density both would be the 
same, but by hypothesis the system lasted so long in state si 
that is at "equilibrium", and only the coarse-grained den- 
sity can any longer be relevant to the dynamics of the fi- 
nal relaxation. We show that the outcomes predicted with 
energy E2 and volume functions 70(77) and Ti(rf) are cer- 
tainly differerrtJn_yieIjynden-Bell theory of violent relax- 
ation llLvnden- Bell 1967). as well as in a more recent theory 
by Nakamura iNakamura 1120001) . 

This paper is organised as follows: in Sec. [5] we give a 
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short overview of some of the difficulties in the theory of 
violent relaxation, which we feel complement the main sub- 
ject of this paper. Then in Sec. [2] we demonstrate the non- 
transitivity of the Lynden-Bell theory of violent relaxation, 
using the thought experiment that was presented above. In 
Sec.^Jwe give a brief description of Nakamura's theory which 
is based on the information-theory approach. We re-derive 
his theory using a combinatorial approach that enables us 
to compare it to Lynden-Bell's theory. Then in Sec. we 
demonstrate that also Nakamura's theory is non-transitive, 
using the thought experiment once again. In Sec. |S| wc 
present our conclusions. 



2 AN OVERVIEW OF THE DIFFICULTIES IN 
THE THEORY OF VIOLENT RELAXATION 

In this section we offer a short discussion of other problems 
connected with the process of violent-relaxation and 'theo- 
ries' that aim to predict its outcome. Let us start from first 
principles. Most large iV-body systems governed by long- 
range forces which are not initially in balance will oscillate 
with decreasing amplitude before they settle into a state in 
which the potential of the long range force becomes almost 
steady. Such violent-relaxation processes are known to oc- 
cur in gravitational iV-body systems. Thereafter evolution 
may continue due to the shorter range interaction in which 
the graininess of the individual particles is of importance, 
but there is a large class of systems in which this secondary 
evolution is on a much longer timescale. Violent relaxation 
under gravity does not last long. After a few oscillations on 
the timescale (Gp) _1,/2 it is over. Thus the whole idea that 
the interaction of the particles with the mean field will lead 
to some unique detailed statistical equilibrium state, only 
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dependent on the initial conditions via the dynamically con- 
served quantities, is more of a vain hope and a co nfession 
of ign orance than an established fact. Nevertheless. lfl6non I 
( 1968) gave some evidence in fa vour of its prediction and 
for co s mological in itial conditions iNavarro. Frenk fc White I 
il995l fl"996|. |l997 l) show a considerable universality in their 
results. iBinnev I (120041) gives a lovely toy model which he 
uses to criticise deductio ns from JV-body sim ulations. Even 
in the initial discussion jLvnden-Bell I 1967) it was admit- 
ted that there were many stable steady states into which a 
gravitating system could settle and that violent relaxation 
would be incomplete so that the system would not inevitably 
attain a state close to the more probable one. 



A second worrying aspect of any equilibrium theory 
is the apparent lack of an analogue to the law of detailed 
balance. At real thermodynamic equilibrium there are no 
cyclic processes going around and around, but each indi- 
vidual emission process is exactly balanced by the corre- 
sponding absorption process. In radiation theory Einstein 
introduced his stimulated emission process just to ensure 
that this would be so. In the process of violent relaxation 
each element is interacting with the potential of the whole 
system, and one might expect some to be highly accelerated 
as Fermi argued to get his cosmic ray acceleration process. 
Conservation of energy must however lead to some dynam- 
ical friction term but this cannot be mass related as that 
is at odds with our earlier arguments that energy gain is 
independent of mass. 



3 NON-TRANSITIVITY IN THE 

LYNDEN-BELL THEORY OF VIOLENT 
RELAXATION 

The Lynden-Bell theory of violent relaxation (hereafter 
LB67) aims at predicting the final equilibrium state of a 
collisionless gravitating system undergoing a violent relax- 
ation. In order to demonstrate its non-transitive nature we 
shall first recall its main results for the general case where 
initially the system has more than one phase-space density 
levels. 



3.1 Main results of the LB67 theory 

In the LB67 theory we deal with a system which is initially 
out of equilibrium. Its initial state is specified by its to- 
tal energy E and the phase-space volumes Vi, V2,. . . of the 
initial phase-space density levels 771, r)2, - • - 1 . Phase-space is 
then divided into micro-cells of fixed volume uj, which can 
be either empty or hold a phase-space element of one of the 
prescribed levels. The state of all these micro-cells defines a 
micro-state. 

Next, we let the system re-distribute its phase-space 
elements as it approaches an equilibrium. The micro-cells 
are then grouped into macro-cells, each macro-cell contain- 
ing v micro-cells. A macro-state of the system is defined 
by the matrix {riij} which specifies how many phase-space 
elements of type J ended up in the macro-cell i. The coarse- 
grained phase-space density function (DF) at the macro-cell 
i is therefore given by 



Si 



1/ — ^ 



(1) 



Finally, not all systems have violent relaxation. In 
some earl y experiments with pulsating concentric spheri- 
cal shells, iHenon I jl968f) found a few examples of grav- 
itating system with persistent oscillation that defied the 
general decay. Newton in Principia showed that systems 
with a force law between particles proportional to sepa- 
ration (rather than inverse square) oscillated forever. In- 
deed Newton solved that iV-body problem completely, show- 
ing that each particle moved on a central ellipse centred 
on the barycentre and that all those orbits had the same 
period. Lynden-Bell & Lynden-Bell (1999a) showed that 
Newton's work could be extended to all systems in which 
the total potential energy was any function V(r) where 
r 2 = EiM^i ~ rf/M = lEi^^-fc - ^f/M and 
r is the position of the barycentre. Likewise, the iV-particle 
Schrodinger equation was exactly sol ved for such potentials 
l)Lvnden-Bell fc Lvnden-Bell Ill999b^ : the most interesting 
example of which is probably V = GM 2 /r. This has a 
force between any two particles of Gm,imj(ri — r,)/r 3 with 
the final r 3 not \n — | 3 but the r defined above. A re- 
lated exception to violent relaxation with a most interesting 
ever-pulsating detailed "Maxwellian" statistical equilibrium 
is the system with an inter-particle potential iwiiTJij (jto"y + 
k'r~ 2 ). In spite of the complication of the individual particle 
motions due to the inverse cubic repulsion between particles, 
such systems pulsate i n scale forever in exact undamped si m- 
ple harmonic motion dLvnden-Bell fc Lvnden-Bell 1 12004 ) . 



Then, in the spirit of ordinary statistical mechanics, one 
assumes that the system has an equal a priori probability of 
being in each one of the micro-states. To find the equilibrium 
state one maximises the function W^({nij}) which counts 
the number of micro-states that correspond to the macro- 
state {riij}, hence obtaining the most probable state. 

When maximising W one has to consider only those 
macro-states for which the total energy is E and the overall 
volume in each of the initial phase-space levels is Vi, V2, ■ ■ ■■ 
This can be done in the usual way with Lagrange multipliers. 
After passing to a continuous description of the macro-cells, 
the resultant DF is 

/Ls(r,v) = / ifl [e(r,v)] = A(e) £ „ jB -*w(«-w) , (2) 

j 

with 



A(e) 



1 



(3) 



1 + Y^je-^J^-^j) 

Here e(r, v) = v 2 /2 + $Lfl(r) is the energy per unit mass of 
the (r, v) phase-space cell, <&lb(t) is the gravitational po- 
tential, calculated self-consistently from the Poisson equa- 
tion 



V 2 $(r) = 4ttG Jd 3 vf LB (r,v) 



(4) 



1 In the original formulation of LB67, Lynden-Bell used the 
masses of the different levels instead of their phase-space volume. 
However these are trivially related to each other by Mj = rjjVj 



Inconsistency in theories of violent-relaxation 3 



The dimensional constants /3, fii, [12, . . . are the Lagrange 
multipliers, which are calculated from the energy conserva- 
tion constraint 



Jd\f LB (r,v) 
itial conditions 

Jd\A{e)e- l3r]A '-' lj) =Vj, J 



2 2 

and the initial conditions constraints 



= E 



= 1,2, 



(•>) 



(6) 



where we have used used the notation d r = d rd v. 

Equations 12161 arc the main results of the LB67 theory. 
For our needs, however, two small modifications are needed. 
Firstly, by introducing the energy density function 



dV5[e(r,v)-e ] 



(7) 



the integral in the LHS of Eq. Jfjjl can be written as 
J<?(e)A(e)e~' 3 '' J ' E ~'' J ' de. Secondly, we pass to a continu- 
ous description of the initial density levels using the phase- 
space volume distribution function T(r/), which is defined so 
that T(r))dr) is the phase-space volume initially occupied by 
phase-space densities in the range 77 — > r\ + drj. Formally, if 
fi (r, v) is the initial DF then t(t]) is given by 



t(ij)= /d 6 T<5[/,(r,v)-,7] 



(8) 



By letting each density level rjj have a small width A 77, the 
j sums in Eqs. |3J can be changed to integrals by 



E 



hence Eqs. ||2T51 are now 
1 



A(e) 
V 2 $(r) 
E 



A77 



A(e) 



Ar) 



-0V(t-Hri) 



drj , 



drj 



AnG /d 3 v/ LS (r,v) 



v 1 / % 



d T/ LB (r,v 
Jg{e)A{e)e-^ (l -^ ) de 



(9) 

(10) 
(11) 
(12) 
(13) 
(14) 



Notice that the amplitude of Ar] is unimportant as it can 
always be absorbed into the Lagrange multipliers jj, v . 

Finally we recall that in a spatially infinite domain 
equations 11011411 have no solution since the density can 
spread indefinitely while conserving its energy and increas- 
ing its entropy. A common way to overcome this problem, 
which shall be adopted here, is to work within an rigid sphere 
of radius R, and to assume that the resulting equilibrium 
configuration is spherical. Such model, although not realis- 
tic, is an easy way to obtain a finite solution. 



3.2 The double relaxation experiment 

To test the transitivity of the LB67 theory we propose the 
following four-steps thought experiment which was men- 
tioned in the introduction: 



(i) We prepare a system with one density-level (the water- 
bag configuration) and a total energy of E. In accordance 
with the above section, we put the system inside a sphere of 
radius R. We denote the initial state of the system by so. 

(ii) We let the system go trough a violent relaxation pro- 
cess to an equilibrium which is denoted by s\. 

(iii) We add an amount of AE energy to the system by, 
for example, a strong impulse of an external gravitational 
field. The system then goes once again through a violent 
relaxation process and settles down in a new equilibrium 
state, denoted by S2. 

(iv) We prepare a new system with the same parameters 
as so except for the energy, which is set to E + AE. We let 
it go through a violent relaxation process to an equilibrium 
which is denoted by S3. 

Now if the theory is transitive then necessarily S2 = S3. 

To see if this is really the case in the LB67 theory, 
we begin by calculating the si and S3 states which are 
relatively simple to calculate, being the outcome a water- 
bag configuration. This calcu lation has been fully done in 
IChavanis fc Sommeria I Jj.998), and here we use some of their 
results. 

The so step was prepared with total mass M = 1, G = 
1, R = 1 and an initial phase-space density level 

10 3 

770 = ; , (15) 

V5l2iv 4 G s MR 3 

which according to IChavanis fc Sommeria I <ll998f) guaran- 
tees that for each energy there would be only one equilib- 
rium state. 

From Eq. © we find that the DF of si and S3 is the 
well-known Fermi-Dirac distribution 



AW = 
/ 3 (e) = 



rjo 



I -|- e /3i?7o(E-Ml) 



J/0 



1 + e^oO-M) 



(16) 
(17) 



with j3i,j3z, /Ji, /J3 Lagrange multipliers to be fixed from the 
energy constraint and the initial conditions. As there is only 
one density- level, the initial condition constraint can be re- 
placed with the conservation of mass constraint. 

In IChavanis fc Sommeria I il998l) it is shown how the 
Lagrange multipliers can be found for a given mass and en- 
ergy, and we therefore do not repeat these steps here but 
instead give the values of these parameters in Table 1. All 
numerical calculations were done using the GNU Scientific 
Library 1.5 (GSL 1.5), which is a free software available 
from http://www.gnu.org/software/gsl/. The differential 
equation for <J?(r) was solved using an embedded Runge- 
Kutta Prince-Dormand (8,9) method, whereas integration 
was done using a 51 points Gauss-Kronrod rule. In all cal- 
culations a relative error of less than 10 -3 was maintained. 

Figure shows the radial density profiles of si and S3. 
As the S3 configuration is more energetic, its mean kinetic 
energy is higher and as a result the distribution is less con- 
centrated than the si distribution, and has a lower density 
core. Figure shows the DFs of si and S3. Notice how both 
distributions have a substantial degenerate part, as for both 
configurations the Fermi energies are e ai ~ $i(0)/2. 

Finally, to verify that the S3 configuration is indeed 
more mixed than the Si configuration - and therefore a 
transition si — > S3 is permitted by the mixing theorem 
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Figure 1. Density profiles of the si and S3 states in the double- 
relaxation experiment of the LB67 theory, p is the average density 
given by p = 3M/ (4nR 3 ). The S3 state corresponds to a hotter 
system with higher energy which makes it less concentrated. 
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Figure 2. The phase-space densities of si and S3 as functions of 
the energy per unit mass e. The functional form of these densities 
is the well-known Fcrmi-Dirac distribution, given by Eqs. I16II17I 
and Table. 1 



jTremaine. Henon fc Lvnden-Belll Il98fl - we have calcu- 
lated the M(V) functions of Si and S3. The M(V) function is 
defined in the following way: we first define M(rj) and V(rj) 
as the cumulative mass and phase-space volumes above the 
phase-space density 77, 



M(n) 



d\ H[f(r,v) - r,] = 



r<R 



v' T (v') drf 



(18) 
(19) 



with H(-) being the Heaviside step function. As V(rj) is a 
monotonically decreasing function, we can invert it and ex- 
press r? in terms of the cumulative volume V . Plugging this 
into M(rj) we get the M(V) function. According to the mix- 
ing theorem, the distribution S3 is more mixed than the 
distribution si if and only if Ms(V) ^ Mi(V) for every V, 
as is shown in Fig. [3] 



3.3 Analysing the S2 configuration 

Let us now turn our attention to the S2 configuration. Seem- 
ingly, we need to calculate the n(rj) function of the si con- 
figuration, and together with an energy of E + AE, solve 




Figure 3. The M(V) functions of si and S3. Mo = 1 is the 
total mass and Vq = 1 is the initial phase-space volume that is 
occupied by the density level rjo. We notice that M^IV) $C M\(V) 
for every V and therefore S3 is more mixed than si according to 
the mixing-theorem. 



the equations 
Me) = ^A(e) 
A(e) : 



1 f 00 

Ar l J*,(0) 



(20) 
(21) 

(22) 



The function 32(e) needs to be calculated from Eq. JjJ using 
the gravitational potential $2(7"), which has to be recovered 
from /2(e) using the Poisson equation 111 2t . Finally, the re- 
sultant /2(e) would be compared to /3(e) to see if the two 
configurations are equal. 

There is, however, a much simpler way to see if /2(e) = 
/3(e). Let us assume that indeed this is the case, and 
that consequently also $2(7") = $3(r) and 52(e) = 53(e). 
In such case it is possible to recover the full expression 
-^A(e)e~ l32,lo( - e ~' 1 i ) in terms of the known functions /3(e) 
and 53 (e). 

We start by replacing /2(e) «-> /3(e) and 52(e) <-> 53(e) 
in Eqs. 1201221 . Differentiating Eq. 1211 with respect to e 
and substituting /3(e) from Eq. 120H . we obtain 



A 2 (e) = frA(e)Me) , 

which yields 

A(e) = a exp I —fa / / 3 (e')de' I = a Q A (e) 



(23) 



(24) 



Here ao is an unknown integration constant. The integral in 
Ao(e) can be easily done analytically if we recall the defini- 
tion of /3(e) which is given in Eq. 1171 . yielding 



lo(e) = [ 



1 + e 



-031o( £ -M3)l ^3 



(25) 



Having found A(e), we use Eq. i'2'It to find the Lagrange 
multipliers fj, v : 



1 r 00 

[nM] ao^- / 53(e)A (e) C - & " e de. (26) 



Inconsistency in theories of violent-relaxation 5 



Table 1. The numerical parameters that specify the si and S3 states. E is the total energy of the state, 7)0 is the initial density level 
given by Eq. 1151 . /3,fi are the Lagrange multipliers in Eqs. 1161 1171 . and ( I > (0) is the gravitational potential at r = 0. 



State 


E 




/' 


3.(0) 


si 


-6.157 


1.2 


-14.87 


-33.96 


S2 


-1.589 


1.0 


-8.322 


-15.12 



Therefore 

A{e) 



Arj 



Bo{t)) 



Me) 



with 



Bo(v)= / 93(e) A (e')e-^ e ' de' 

J$ 3 (0) 



This finally gives us 

/ 2 (e) = A (e) 



3 71 (77) -force. 



(27) 



(28) 



(29) 



Note that the unknown integration constant 00 has been 
cancelled out. The only remaining unknown is 02 which can 
be fixed by requiring that the energy of S2 will be equal to 
E + AE. Once this is done, we have an expression for /2(e) 
which is equal to /3(e) if and only if S2 is identical to S3. 

The procedure above is mathematically straightfor- 
ward, however, numerically it is slightly more complicated 
as the Ti(rf) has a very strong peak near r\ — rjo due to the 
degeneracy. It is therefore preferable to perform the calcula- 
tion using the cumulative version of 71(77), which was defined 
in 11811 . For a spherical, isotropic system in a sphere of ra- 
dius R with a DF /(e) and a gravitational potential <&(r), it 
is easy to verify that V{rj) is given by 



1/(77) 



,3/2 



(4rr) 2 



s"[e(77)-$(s)] 3/2 ds 



(30) 



with e(rj) being the inverse function of /(e), and r(r7) is 



R 



e(77) < 



-GM/R 
-GM/R 



(31) 



Once Vi (77) is calculated from the formula above [using /1 (e) 
and <3>i(e)], we can calculate /2(e) from Eq. I)29|l using inte- 
gration by parts: 



f 2 (e)=A (e) / ^1(77) 



l-/3 2 7?e- 



(32) 



K3 



GO 



40 



20 




0.5 



1.5 



Figure 4. The energy E2 of different values of /32, together with 
E 3 = -1.589. P2 = 0.3707 is the value that gives E 2 = E 3 . 
For higher values of up to P2 = 100 we did not find any other 
solution and therefore we believe that the above solution is the 
only physical solution. 
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Figure 5. Comparing /3(e) to /2(e) which was derived using 
Eq. 1321 . The two functions are unequal proving that the S3 state 
is different from the S2 state. 



3.4 Results 

To satisfy the energy constraint E2 — E3 — —1.589, we 
calculated E2 for various values of fa and chose 02 = 0.37071 
which gives the correct energy as shown in Fig. [1] We did 
not find any other solution in the range 0.1 < (3 < 100 and 
therefore we believe that p2 = 0.37071 is the only relevant 
solution. 

Figure shows the graphs of /2(e) and /3(e) once 02 
was fixed. Clearly the two graphs strongly disagree, in some 
places by more than one order of magnitude - much more 
than the numerical error in our calculations. The conclusion 
is therefore that the LB67 theory is not transitive. 



4 THE INFORMATION-THEORY APPROACH 
TO VIOLENT RELAXATION AND ITS 
RELATION TO THE LB67 THEORY 

Recently, a new approach to violent relaxatio n was pro- 
posed in a interesting paper by T. K. Nakamura dNakamura I 
2000). In that paper, Nak amura uses an information-theory 
approach llJavnes Hl957aHbh to define the entropy of a colli- 
sionless system and thereby find its equilibrium state. Naka- 
mura's theory (hereafter NK00) predicts a different equilib- 
rium state than LB67, and it is therefore interesting to check 
weather his theory is transitive or not. We will not, however, 
try to answer the question which one of these theories is 
more correct as it is, in our opinion, still an open question. 
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Instead, we shall first give a brief description of NKOO and 
its main results, and then re-derive it theory using a com- 
binatorial approach which would enable us to compare it 
with the LB67 theory, and point to the reasons of why they 
differ. Finally we will analyse the two-levels configuration 
which goes into the water-bag configuration in a limiting 
case. The result of this analysis will be used in the next sec- 
tion when we examine the transitivity of the NKOO theory 
in a double relaxation experiment. 



4.1 An outline of the NKOO theory 

In the NKOO theory we adopt the probabilistic description 
of the phase-space density f(t,r,v). Let /o(r, v) be the ini- 
tial phase-space density of the system, and define the initial 
probability distribution po(r,v) for finding a (single) test 
point at t = by 



Po(r,v) 



M 



/o(r,v) 



(33) 



Then we let the test point move under gravity just like any 
phase-space element, and we define the probability distribu- 
tion p(r, v, t) as the probability distribution of finding the 
test point at time t > 0. The conservation of phase-space 
volume guarantees that p(r, v, i) — /(r, v, i)/M for all t > 0. 

Next, we divide phase-space into macro-cells i = 
1, 2, 3, . . . of volume uj, and define the coarse-grained prob- 
ability pi as the probability of finding the point in the i'th 
macro-cell when the system reaches an equilibrium. From 
the above discussion it is clear that pi is equal to fi/M with 
fi being the coarse-grained DF in the macro-cell i at equi- 
librium. 

To calculate pi using the information-theory approach, 
we define the joint probability- distribution pi(r,v) which 
measures the probability of initially finding the test point 
at the (r, v) and later at the macro-cell i. Then we max- 
imise the Shanon entropy 



S=~^2 / rf6r Pi( r > V ) lo gPi( r > V ) 



(34) 



subject to constraints of energy conservation, phase-space 
volume conservation and initial conditions. The resultant 
distribution can be best written in terms of the conditional 
probability Ki(r,v) = Pi{r, v)/po(r, v): 



Ki(r,v) = e -^i-*( r .v)-Ai/ P0 (r,v) 



(35) 



with €i being the energy per unit mass of the macro-cell i, 
and 0, <5(x, v), Xi are the Lagrange multipliers, to be found 
from the energy conservation constraint, from the initial con- 
ditions 

X)tf 4 (r,v) = l, (36) 

i 

and from the phase-space volume conservation constraint 



Ki(v, v) d r = Q 



(37) 



From Eqs. 1351371 it is evident that the dependence of 
Ki(r, v) on the indices i, r, v is only via 6; and po(r, v), the 
latter can be trivially be replaced by /o(r, v). Therefore, the 
above equations can be re- written using the K(e, rj) function, 



together with the t(t]) and g(e) functions which were defined 
in Sec.EHl 

K(e,n) = ^y^ e -^-«M-He)/v ; (38 ) 



r- do 

/ 9 

J 0(0) 



g(e)K(e,r,)de = 1, 



r(ri)K(e, rf) dr\ 



The coarse-grained equilibrium DF is then given by 



(39) 



(40) 



(41) 



As noted by Nakamura, a prominent difference between 
his result and the LB67 results is that in the non-degenerate 
limit his expression reduces to a single Maxwellian distri- 
bution, whereas the LB67 expression is a superposition of 
Maxwellian distributions with different dispersions. This dif- 
ference can be attributed to the fact that in LB67 we discre- 
tise phase-space using phase-space elements of equal volume 
and different masses, while, as we shall see below, the NKOO 
theory can be derived by discretising phase-space using el- 
ements of equal mass, which are associated with different 
phase-space volumes. 

Another evident difference comes from the phase-space 
volume conservation constraint Eq. 1361 . This constraint 
guarantees that the total phase-space volume of all phase- 
space patches that ended up in macro-cell i will be equal to 
the macro-cell volume. Consequently the total phase-space 
volume of the initial system J^°T(r})dri must be equal to 
the total phase-space volume of the non-vanishing phase- 
space density in the equilibrium configuration. This con- 
straint does not exist in the LB67 theory where the macro- 
cells can be only partly full - as is the case, for example, in 
the non-degenerate equilibrium of a system which is initially 
in the water-bag configuration. Furthermore, this can not 
be trivially changed by adding a volume of zero phase-space 
density to the initial condition, because setting p(r, v) = 
would lead to divergences in Eq. 1351 . In Sec. 14.31 we shall 
see how this problem can be overcome by using a limiting 
procedure. 

4.2 Deriving the NKOO theory in a combinatorial 
approach 

To derive the NKOO theory in a combinatorial approach, we 
realise the phase-space density distribution using N 3> 1 
elements of equal mass m. As in Sec. 13.11 we assume that 
initially the system is made of a discrete set of density levels 
7/1,772,... occupying phase-space volumes Vx,Vi,.... Then 
the overall number of elements that realise a phase-space 
density nj is Nj — Vjnj/m. 

Next, we let the system reach an equilibrium through 
the process of violent relaxation, and divide phase-space into 
macro-cells of equal volume Co, which are label by the in- 
dex i = 1, 2, We define a micro-state by specifying the 

macro-cell in which every element ended up. A macro-state 
is then defined by the matrix {nij} which counts how many 
elements that initially realised the density level r\j ended up 
in the macro-cell i. Using {nij}, the coarse-grained DF at 
macro-cell i is given by 

fi (42) 
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Finally, we define the function W[{nij"j which counts 
how many micro-states give the macro state {riij}. It is then 
a simple combinatorics to show that 



W({nu})=H 



U,"ul 



(43) 



Let us pause here and explain that this rather simple for- 
mula is a result of the way we define a micro-state - by 
specifying the macro-cell in which each element is found. 
We do not care exactly how the different elements are dis- 
tributed in each macro-cell. A different approach, in which 
the macro-cells can be only partly full, and the distribution 
of the different elements in a macro-cell is taken into accoun t 
when defining a micro-state, was taken bv lKull et al I il997f) . 
It is not difficult to see that when one adds the constraint 
that all macro-cells must be completely full to their theory, 
Nakamura's results are recovered. This is because in such 
case the number of different ways to arrange the different 
elements in the macro-cell is independent of which elements 
we are organising - as long as the macro-cell is completely 
full. Therefore the number of micro-states in a macro-state 
would be proportional to Nakamura's W({riij}), and con- 
sequently the equilibrium state would be identical in both 
theories. 

Next, we use W({riij}) to define the entropy 



S = log W ~ const — S~] n ij(}°e n iJ — 1) 



(44) 



where in the second equality we have used Stirling's formula 
to approximate log(nij!) ~ raj (log n;j — I). 

Before maximising the entropy to find the most prob- 
able macro-state, we first write down the constraints on 
{riij}. The first constraint comes from the initial conditions 



Vjrjj 



(45) 



Then we have the phase-space volume conservation con- 
straint, ensuring that the total phase-space volume that is 
carried by elements that ended up in the macro-cell i will be 
exactly u), or, in other words, that each macro-cell is com- 
pletely filled: 

at — > m 

> rnj — = oj . (46) 
j VJ 



The last constraint is the energy constraint 



E^ 



= E 



(47) 



where n and Vi are the mean position and velocity of the 
i'th macro-cell. 

To maximise the entropy under the above constraint 
we use Lagrange multipliers. The function that we wish to 
maximise with respect to n^j is therefore 



1 = 8— ^2 ^i n iJ — E Sjntj — (3E 



(48) 



Differentiating / with respect to rnj and equating it to zero, 
we get 

dl Xi 



driij 



log nu 



VJ 



f3mei — 



(49) 



with Si = vf /2 + <$>(ri) as usual, and therefore 



— 3me,- — 5 j — 

riij = e J "j 



(50) 



Finally, we pass to a continuous description by giving 
every initial phase-space density level a small width An. 
Then using the r(r]) and g(e) function we replace 



E 

j 

E 

i 

Vj 
nu 



I 

An 

1 



dn 



deg(e) 



ArjT(nj 
Anil) 



T(nj)njK(ei,nj) 



(51) 

(52) 

(53) 
(54) 



Plugging these replacements into Eqs. I5UI 1451 14611 and re- 
defining m/3 -> /3, we recover the NK00 Eqs. QHlEHEOt- 

This combinatorial formulation of the NK00 theory is 
very much along the lines of ordinary statistical mechanic 
of a classical Boltzmann gas. Indeed, if we replace the no- 
tion of phase-elements with particles of equal mass and dis- 
card the constraint of conservation of phase-space volume 
Eq. 1401 . we have a text-book derivation of the Boltzmann 
gas statistics. It is therefore not surprising that Nakamura 
found that his equilibrium DF reduces to the well-known 
Maxwell-Boltzmann distribution in such case. We do not 
agree, however, with Nakamura's claim that this property 
is a proof for its correctness over the LB67 theory. This is 
because a collisionless relaxation is essentially a very differ- 
ent process from the collision-full relaxation that occurs in 
Boltzmann gas, driven by different physical processes over 
different timescales. However, as previously mentioned, de- 
ciding which theory is more correct is not the goal of this 
paper. 

4.3 Analysing the two-levels configuration 

As was noted in the end of Sec. 14. II the NK00 theory cannot 
handle a zero phase-space density directly. Therefore it is not 
straightforward to analyse the equilibrium state that results 
from an initial water-bag configuration, as in this configura- 
tion there is one patch of phase-space density 770 surrounded 
by an infinite volume of zero phase-space density. The way 
this can be done is to consider an initial state with two den- 
sity levels 770 and ni with corresponding volumes Vb and V\ . 
The water-bag configuration is then recovered by taking the 
limit rji — > 0, Vi — » 00 and r/iVi — > 0. 

To derive the equilibrium configuration of the two-levels 
system we use the fact that in this particular case the ma- 
trix {nij} can be expressed in terms of fi, thereby greatly 
simplifying the end result. Let us then re-derive the equi- 
librium equation for this particular case instead of using 
Eqs. 138I4H . Denoting by n^o and m,i the total number of 
elements of rjo and rji that end-up in the i'th macro-cell, the 
coarse-grained DF is given by 



fi 



(55) 



Then using the conservation of phase-space volume con- 
straint 1461 1 . 



ni.o 



Nn 



+ UiA 



(56) 
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with No,Ni being the total number of elements with den- 
sities T70 , 771 , given by Eq. 145H . we express n^o and rc^i in 
terms of fv. 



m vo - Vi 



Vo 



& Vo — fi 
m,i = r)\ 



(57) 



(58) 



m rjo - 771 

The energy constraint is given by Eq. 14711 . and the initial- 
condition constraint is 

»7t>Vb 



(59) 



Notice that we need only the Nq constraint since the Ni con- 
straints follows directly from requiring that the total phase- 
space volume occupied by the equilibrium system would be 
equal to V — Vo + Vi. In fact, instead of Eq. 159H . we can use 
an alternative total mass constraint, provided that the over- 
all phase-space volume is conserved. This is done as follows: 
expressing n,,o in terms of fi in Eq. 1591 we get 



Vo - 7]1 

which gives us 

^ ufi — rji ^2 & = V0V0 — ViTji 



(60) 



(61) 



But & = Vo + Vi (conservation of total phase-space vol- 
ume) and therefore we find 



ufi = Vorjo + Vir)i = M 



(62) 



Adding these constraints together with the appropriate 
Lagrange multipliers to the entropy, the expression that we 
need to maximise is 

I ~ const — ^2 n i,o ( log riifi — l) (63) 

i 

i i 

Differentiating / with respect to fi and equating to we 
find 



Vo w ( & fi-r]i 

■ log Vo 



Vo - Vi m 



m rjo — 771 



(65) 



Vi w / ' Cb rjo - fi . 
H log [ — —771 + /3wei + /Hoi = .(66) 



770 - 771 m 



m 770 - 771 



After a trivial algebra and redefinition of the Lagrange mul- 
tipliers /3 and jl to fj and /j, we obtain 



fi-m 
vo-vi 



VO-fi 



m /vo 



(67) 



v >jo -11 , 

Notice how the denominator provides an upper cut-off for 
fi, as it forbids it from exceeding 770. 

Consider now the 771 — > limit. Seemingly, it would go 
into an isothermal sphere 

fi = Vo e- p ^-^ , (68) 

but this is not the case as for every finite 771 , fi cannot exceed 
rjo- It is easy to see that the right limit is therefore 



e < n 
e > n 



(69) 



This distribution is not the LB67 Fermi-Dirac distribution 
given by Eq. I16H or Eq. I17H . but is what corresponds to 
that distribution on the NK00 theory. It is not smooth, and 
is exactly isothermal for energies e > fi. 

For the water-bag model in LB67 the condition that no 
two elements of phase-density can overlap leads to a statis- 
tics with exclusion, equivalent to the Fermi-Dirac problem. 
It is not clear to us how Nakamura's formulation could ob- 
tain the Fermi-Dirac statistics. 



5 NON-TRANSITIVITY IN THE NAKAMURA 
THEORY OF VIOLENT RELAXATION 

Having found the equilibrium configuration of the water- 
bag initial configuration in the NK00 theory, we are in a 
position to test the theory's transitivity. The procedure for 
that is identical to the one that was used in the LB67 case, in 
sections 13. 2113. 31 and therefore will not be repeated. Instead, 
we shall first describe how the si and S3 configurations are 
found and then how S2 is compared to the S3 configuration. 

To find the Si and S3 configurations, we must first find 
the gravitational potential of the DF in Eq. (16911 in a sphere 
of radius R, and then fix f3 and /_t so that the overall energy 
and mass will be equal to E and M. Additionally, just as in 
the LB67 case, we assume that the final equilibrium state is 
spherical and therefore the Poisson equation for $(r) is: 

hi { r2 t) = i**f<*j~*S»*[*l* + *<r)] dv . 



Passing from $(r) to the dimensionless 7/>(r) by 

V>(r) = [3[n- *(r)] , 

Eq. I7(ty simplifies to 

Id/ 2 d^\ , . 2 
-^d-r{ r drj= (47r) GV ° 

, f3v 2 /2 < V»(r) 



(70) 



(71) 



(72) 



X /o v2dV \ e*(V^ i 7v*/2^U (73) 



Finally, changing variables r — > x 

I67TV2G770 x 1/2 

/31/2 



(74) 



and using the Error-function erf(x) = -J= f Q x e t dt, the or- 
dinary differential equation for ip(x) is 



-^te( x d*- r m 



(75) 



with 



7(7/0 



2 c 



ip< 
ii> > 



(76) 



To integrate this equation we must first set its initial condi- 
tion. We let ip(0) = ipo be a free parameter, and ip'(0) = 
since in a spherical system the gravitational force vanishes 
in the centre. Then once ip(x) is (numerically) found, we fix 
/1 by requiring that $(i?) = —GM/R. This way we can find 
the gravitational potential, and thereafter the total mass 
and energy for any given (3 and tpo- The last step is to find 
the right (3 and ipo that would give us M and E. 
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Practically, instead of looking for (3, ipo for a specific 
E and M choice in si and S3, we have picked two different 
values of /3 and and two corresponding values of ipo to satisfy 
the total mass constraint. Once ipo is fixed it also defines a 
total energy. We chose f5\ > ftz in order to obtain Ei < Ez- 

Let us now see how the S2 «-» S3 comparison can be 
done. According to Sec. 21 /2(e) is determined by the follow- 
ing set of equations: 

-/3 2 e-S( V )-X(e)/n 



10 



(a>)- 



1 = 



/2(e) = 



gx{e)K(e,rf) de 

*2(0) 

Ti(r])K(e,r))dr) , 

no 

7"i (ri)riK(e, 77) dr\ 





Additionally, we know that /3(e) is given by 



/ 3 (e) 



?7o 
770 e 



-/33( e -/ J 3) 



e < 



(77) 
(78) 

(79) 
(80) 

(81) 



and 33(e), $3(7-) have been found as described above. Assum- 
ing that S2 = S3, we replace these functions with /2(e), 52(e) 
and <£ 2 (e) in Eqs. 178I80H . and differentiate Eq. (JSDJ with 
respect to e. Using Eq. 179H we get 



ys(e) = -A/8(e)-V(e) , 



(82) 



and therefore 

\{e) = C + j3 2 l f 3 ( e ')de' - f 3 (e) = C + X (e) , (83) 



with C some unknown integration constant. The integral in 
Ao(e) can be done analytically, yielding 



Ao(e) = 



[||-l+/32(M3-e)] , e< M3 



Then from Eq. we find that 

/" OO 

J $ 3 (o) 



-/3 2 e-Ao(e)/77 



rfc 



and therefore 

jr( e ,»,) = (s,)-\ 



-l3 2C -S( v )-X(e)/ v 



-/3 2 e-A (E)/J7 



(84) 



(85) 



(86) 
(87) 



Notice how the unknown integration constant C and the 
dimensional constant u> are cancelled out. 

Next we calculate /2(e) using Eq. 18011 and fix 02 such 
that the energy of the system is E + AE. Once /3 2 is fixed, 
K(e, rf) is completely resolved in terms of si and S3 func- 
tions, and we can check if it solves the maximum-entropy 
equations by plugging it into Eq. 17911 . 

Finally, we should note that as in the LB67 case, the 
ti (77) function has a strong peak at r\ — rjo due to the degen- 
eracy. Here, however, this peak is proportional to 5(r) — r]o) as 
/i(e) = r/o for every $i(0) < e < fj,i. The prefactor in front 
of this delta function is Vd eg - the volume of phase-space for 
which $i(0) < e < fix, which can be easily calculated from 
Eq. 1301 . Therefore, to preform the integration over Ti(rj) 
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Figure 6. Density profiles of the si and S3 states in the double- 
relaxation experiment of the NK00 theory, p is the average density 
given by p = 3M/(47r_R 3 ). As in Fig. the S3 state corresponds 
to a hotter system with higher energy which makes it less con- 
centrated. 




Figure 7. DFs of s± and S3 in the NK00 experiment. The func- 
tional form of these DFs is given in Eq. IHHJ and Table. 2. Unlike 
the Fermi-Dirac DFs of the LB67 theory given in Fig. |3] these 
DFs have a sharp transition between a completely degenerate 
core with f = r]o for e < p to an isothermal envelope for e p. 



in Eqs. 1791 1801 numerically, we first calculate the smooth 
contribution which comes from the Ti{rf) with 77 < rjo, and 
then add the delta-function contribution by evaluating the 
integrands at rj = 770 and multiplying them by Vd eg - 



5.1 Numerical results 

As in the LB67 case, the so state was constructed as a water- 
bag configuration with 77 given by Eq. 1151 and M — 1, 
R = 1, G = 1. The si and S3 configurations were then 
chosen as described in the previous sub-section, by fixing (3 
and varying ipo until the total mass constraint was satisfied. 
The main numerical parameters of these configurations are 
summarised in Table 2. 

Figure |S| shows the density profiles of the Si and S3 
configuration. As expected, the hotter system, S3, has a core 
with a lower density than the si system. Figure |7| shows 
the DF of these systems, and Fig. |H| compares the M(V) 
functions of these two states, showing that Ms{V) ^ Mi(V) 
for every V - and therefore the transition si — * S3 is allowed. 

Finally, to compare the s 2 configuration to the S3 con- 
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Table 2. The numerical parameters that specify the si and S3 states in NKOO double relaxation experiment. E is the total energy, /3 
and are the Lagrange multipliers in Eq. 1691 . and $(0) is the gravitational potential at r = 0. 



State 


E 


P 




0(0) 


si 


-6.405 


3.0 


-15.30 


-35.08 


S3 


-4.306 


0.68 


-14.08 


-27.80 




jf 0.8 - / 

^ " / 

0.6 - ; 

0.4 I ' ■ ' ■ ' ■ 

5 10 15 

Hv/v ) 

Figure 8. The M(V) functions of si and S3 in the double- 
relaxation experiment of the NKOO theory. Mo and Vb are de- 
fined as in Fig. |3 but here the y-axis is M/Mq whereas in Fig. [5] 
is ln(M/Mo). This is because the range over which M/Mq varies 
is much smaller than that of the NKOO case. Nevertheless, as in 
Fig. |3 M:j(V) M\(V) for every V, showing that S3 is more 
mixed than s±. 



figuration we have varied 02 in the range 0.1 < P2 < 100 
until we obtained E 2 = E 3 = -4.306 with /3 2 = 0.6944. This 
was the only solution in that range, and we believe that it 
is the only physical solution in general. The graph of i?2(/9) 
verses £3 in the range [0.1,4] is shown in Fig. [5] 

Once P2 was found, we used the expression of K(e, /i) 
in Eq. 11850 . to calculate the RHS of Eq. I|790 and compare 
it to 1. Figure [TU1 shows this comparison. The disagreement 
between the RHS and 1 is sometimes as high as 3 orders 
of magnitude - much higher than any possible numerical 
error. We therefore conclude that also the NKOO theory is 
non-transitive. 



6 CONCLUSIONS 

In this paper we have demonstrated that the statistic- 
mechanical theories of violent relaxation by Lynden-Bell and 
Nakamura are both non-transitive. This non-transitivity is 
a result of the phase-mixing that occurs when the system re- 
laxes; as the fine-grained phase-space density filaments be- 
come thiner and thiner, the system is better described in 
terms of the coarse-grained phase-space density. Any fur- 
ther relaxation of the system should be therefore considered 
in terms of the coarse-grained phase-space density - which 
as we have seen would yield different results from a pre- 
diction that is based on the initial fine-grained phase-space 
density. This is a worrying aspect of these theories as it is 
easy to imagine a scenario where part of the system mixes, 
then fluctuates, and then mixes once again. The predictions 




_10 I 1 1 1 1 1 1 ■ 1 

12 3 4 



Pa 

Figure 9. The energy E2 of S2 as a function of 0%, together 
with the energy E3 = —4.306. The value of 02 for which the two 
energies agree is 02 = 0.6944. No other solutions were founds 
for /3 2 up to 100, and we therefore believe that this is the only 
physical solution. 



- 




-10 1 1 1 ' 1 1 1 1 1 

-30 -20 -10 10 

e 



Figure 10. Comparing the RHS of Eq. 1401 to its LHS (in natural 
logarithmic scale) . The two sides of the equations are in gross dis- 
agreement, proving that also the NKOO theory is non-transitive. 

of the theory, based on the fine-grained density, will then 
give us a wrong result. 

In some sense we have been breaking into an open door. 
Even without considering the non-transitivity of the theo- 
ries, they are plagued by severe problems. There exist two 
equally plausible ways of discretising phase-space, one with 
equal volume elements and one with equal mass elements, 
which yield two different results. More importantly, the abil- 
ity of the theories to predict the final outcome of a violent- 
relaxation process is very limited. Indeed, as was mentioned 
in Sec.|5J the most important reason for this is that violent- 
relaxation is almost never complete; the fluctuations of the 
gravitational potential die much faster for the system to set- 
tle in the most probable state. 
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Nevertheless, we believe that these difficulties and am- 
biguities in exactly how to do the statistical mechanics of the 
collisionless Boltzmann equation teach us an important les- 
son. The non-transitivity that we have shown is a sign that 
a kinetic description of violent relaxation is probably incom- 
plete, as the equilibrium is dependent on the evolutionary 
path of the system. Instead, what is probably needed is a dy- 
namical approach to the problem. Indeed most of the above 
difficulties are circumvented if instead of aiming to derive a 
universal most probable state, we reduce our aim to that of 
finding an appropriate and useful evolution equation for the 
coarse-grained /. 

An interesting attem pt to find such equation was 
taken bv lChavanis I il998f) . who used the maximal entropy- 
production principle (MEPP) to obtain a close equation for 
/. His analysis, however, uses the initial fine-grained r(rj) 
to define the (Lynden-Bell) entropy rather than the instan- 
taneous, coarse-grained T(rf), which according to the above 
discussion is more correct. Derivation of a useful dynamical 
equation for / thus remains a challenging open problem. 



ACKNOWLEDGMENTS 

We thank Peter Johansson for his help with the manuscript. 
This work was supported by a Marie-Curie Individual Fel- 
lowship of the European Community No. HPMF-CT-2002- 
01997. 



REFERENCES 

Binney J., 2004, MNRAS, 350, 939 
Chavanis P.-H., Sommeria J., 1998, MNRAS, 296, 569 
Chavanis P.-H., 1998, MNRAS, 300, 981 
Jaynes E. T., 1957, Phys. Rev., 106, 620 
Jaynes E. T., 1957, Phys. Rev., 108, 171 
Henon M., 1968, Bull. Astronomique, 3, 241 
Lynden-Bell D., 1967, MNRAS, 136, 101 
Lynden-Bell D., Lynden-Bell R. M., 1999a, Proc. R. Soc. 
A, 455, 475 

Lynden-Bell D., Lynden-Bell R. M., 1999b, Proc. R. Soc. 
A, 455, 3261 

Lynden-Bell D., Lynden-Bell R. M., 2004, J. Stat. Phys, in 

press, (cond-mat/0401093) 
Kull A., Treumann R. A., Bohringer H, 1997, ApJ, 484, 58 
Nakamura T. K. 2000, ApJ, 531, 739 
Navarro J. F., Frenk C. S., White S. D. M, 1995, MNRAS, 

275, 720 

Navarro J. F., Frenk C. S., White S. D. M, 1996, ApJ, 462, 
563 

Navarro J. F., Frenk C. S., White S. D. M, 1997, ApJ, 490, 
493 

Tremaine S., Henon M., Lynden-Bell D., (1986), MNRAS, 
219, 285 



